Quantum Critical Behavior in Kondo Systems 
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This article briefly reviews three topics related to the quantum critical behavior of certain heavy- 
fermion systems. First, we summarize an extended dynamical mean-field theory for the Kondo 
lattice, which treats on an equal footing the quantum fluctuations associated with the Kondo and 
RKKY couplings. The dynamical mean-fleld equations describe an effective Kondo impurity model 
with an additional coupling to vector bosons. Two types of quantum phase transition appear to 
be possible within this approach — the flrst a conventional spin-density-wave transition, the second 
driven by local physics. For the second type of transition to be realized, the effective impurity model 
must have a quantum critical point exhibiting an anomalous local spin susceptibility. In the second 
part of the paper, such a critical point is shown to occur in two variants of the Kondo impurity 
problem. Finally, we propose an operational test for the existence of quantum critical behavior 
driven by local physics. Neutron scattering results suggest that CeCus-xAua; passes this test. 
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PACS numbers; 71.10. Hf, 71.27.-Ha, 71.28.-Hd, 74 
I. INTRODUCTION 

Heavy fermions close to a quantum phase transitionB 
represent an important class of non-Fermi liquid metals. 
The transition is typically between a paramagnetic metal 
and a magnetic metal. The generic phase diagram is illus- 
trated in Fig. ^, where S represents some tuning param- 
eter. Non-Fermi liquid behavior occurs in the quantum 
critical regime-about 6 = 6c- In one class of materials, 
which includesfi CePd2Si2 and Celna, increasing S corre- 
sponds to decreasing external pressure. In other cases, S 
controls the stoichiometry. For instance,u CeCue-xAux 
is paramagnetic for x < Xc ^ 0.1 and antiferromagnetic 
for X > Xc- 

The magnetic phase transition in heavy-fermion sys- 
tems is generally thought to result from competition be- 
tween Kondo and RKKY physics. The Kondo interac- 
tion tends to quench the local moments using the spins 
of the conduction electrons, whereas the RKKY inter- 
action promotes local-moment ordering. The transition 
occurs when these two processes are of roughly equal im- 
portance. 

Historically, the interplay between Kondo and RKKY 
effects was first studied in the "Kondo necklace," a 
one-dimensional Kondo lattice model simplified by the 
suppression of all charge degrees of freedom. Using a 
(static) mean-field approach, DoniachH found that in- 
creasing /rkky/2k, the ratio of the effective RKKY and 
Kondo scales, leads to a continuous transition from a 
Kondo- insulator-like state to an antiferromagnetically or- 
dered state. Much subsequent work has been |devoted to 
the full Kondo lattice model in one dimension.El A combi- 
nation of numerical and analytical methods have yielded 
results very different from Doniach's: at half-filling, a 
spin liquid forms for arbitrary values of /rkky/Tx; at 
other fillings, the system evolves from a Luttinger liquid 
to a ferromagnetic metal as the ratio Ikkky /Tk is de- 
creased. While these differences stem in part from the 
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special properties of interacting fermions in one dimen- 
sion, they also highlight the importance of quantum fluc- 
tuations in the Kondo lattice. 

The competition between Kondo and RKKY physics 
has alsq-i been studied in the two-impurity Kondo 
problemJj For weak antiferromagnetic RKKY interac- 
tions, single-impurity Kondo behavior is recovered. For a 
very large Ikkky /Tk, by contrast, the two local moments 
lock themselves into a singlet, and no Kondo screening 
takes place. Under special conditions, an unstable non- 
Fermi liquid fixed point separates the two stable Fermi 
liquid regimes. 

In order to better understand the physics of heavy- 
fermion metals, it is necessary to address the dynami- 
cal competition between the Kondo and RKKY interac- 
tions in spatial dimensionalities D n'r-k- In this connec- 
tion, a recently developed extensionl3"ll3 of the dy.n arnical 
mean-field theory (DMFT) of the large-Z? limitEiy of- 
fers a promising avenue of approach. Unlike the standard 
DMFT, which incorporates only purely local dynamics, 
the extended DMFT takes into account nonlocal quan- 
tum fluctuations such as might arise due io RKKY in- 
teractions. The extended DMFT amountallo a resum- 
mation of a diagrammatic 1/D expansion £3lij 

The extended DMFT represents the Kondo lattice by 
a self-consistently determined impurity problem in which 
a single local moment couples to vector bosons as well as 
to a fermionic band. As will be explained below, the 
extended DMFT suggests the possibility of two differ- 
ent types of quantum critical point in the Kondo lattice. 
A scenario of particular interest is a transition driven 
by anomalaus local Kondo physics. We note that other 
authorsE3£3 have also emphasized the role played by local 
dynamics in heavy fermions near quantum criticality. 

The remainder of the paper is organized as follows. 
Section ^ describes the extended DMFT for the Kondo 
lattice. In Section HI, we summarize recent work on 



phase transitions in related Kondo impurity problems. 
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FIG. 1. Generic phase diagram of a heavy-fermion metal 
which exhibits a quantum critical point at temperature T = 
and at a critical value of some tuning parameter, S = Sc. 



Finally, Section IV provides new insights into the quan- 
tum critical behavior in heavy fermions. In particular, 
we introduce a criterion which can be used to identify 
quantum critical behavior driven by local physics. 



II. THE EXTENDED DYNAMICAL 
MEAN-FIELD THEORY 

For simplicity, we neglect valence fluctuations and take 
as our starting point the Kondo lattice model: 

H = ^ Ujcl^Cja- + ^ J Si ■ Sc,i + ^ lijSi ■ Sj, (1) 
{ij),cr i (ij) 

where Si denotes a localized spin on site i, and Sd rep- 
resents the net conduction-electron spin at site i. Equa- 
tion (0) includes both a local Kondo coupling J and an 
explicit exchange coupling between local moments. 
We will consider only the nearest-neighbor terms for lij 
and for the electron hopping tij . (Further- neighbor cou- 
plings can be easily incorporated.) 

Even for = 0, the Kondo coupling in Eq. (Q) gen- 
erates an indirect RKKY interaction between nearest- 



neighbor local moments of strength /: 
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However, this interaction does not survive in rtiie stan- 
dard large-D DMFT for the Kondo latticejid which 
rescales the nearest-neighbor hopping so that — > 

to/^/D, where to is a constant. With this rescahng, 
C!C 1/D, which vanishes faster than the hopping and 
becomes negligible for EL-r^-po. 

The extended DMFTUlil3 ensures that the quantum 
fluctuations associated with the RKKY interaction are 
retained at the mean-field level by rescaling the hopping 
and the explicit exchange coupling with the same power 
of 1/D: t^ij^ to/v^ and /(y) Iq/^/D, where /o is 
a constant. This procedure is well-defined for D ^ oo 
provided that the Hartree contribution is treated with 
care. In terms of diagrammatic perturbation theory, it 
can be shownlj that the extended DMFT is equivalent 



to a conserving resummation to infinite order of a 1/D 
expansion series. 

The mean-field equations resulting from the extended 
DMFT can be expressed in terms of an effective action 
for a single lattice site (the "impurity" site): 

/•/3 rl3 r 

^imp--/ dT dr' ^4(T)Goi(r-r')c.(T') 

Jo JQ [ ^ 

+Sir)-Xo\r-T')S{r') 



+ / dTjS{T)-S,{T)+Si 
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(2) 



Here S'top describes the Berry phase of the impurity 
spin. The Weiss fields Gq ^ and Xo ^ are determined self- 
consistently from the local lattice Green's function Gioc 
and the local lattice spin susceptibility xioc, respectively. 

The effective action of Eq. (g) can be reformulated in 
terms of an impurity Hamiltonian, 
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where the effective dispersion is determined by Gq^, 
while and g are determined by Xo ^- The conduction 
band describes at the single-particle level the effect on 
the selected site of all electrons at other sites in the lat- 
tice. The vector boson fields 4>q extend this description 
to the particle- hole (spin) level by keeping track of nonlo- 
cal quantum fluctuations. These fields are the dynamical 
analogs of the Weiss field in the mean-field treatment of 
classical spin systems. 

Just as in the standard large- 13 limit, once the local 
problem is solved for the impurity self-energy Y.i^p(Lu), 
one can calculate the lattice Green's function via the re- 
lation 



G(k,c.) 



(4) 



where ek describes the dispersion on the physical D- 
dimensional lattice. The lattice spin susceptibility is 
more complicated to determine because it satisfies an 
integral (Bethe-Salpeter) equation rather than an alge- 
braic (Dyson) equation. In the extended DMFT, how- 
ever, the momentum-dependent spin susceptibility is still 
relatively simple: 
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Miio) ~ /(q) ■ 



(5) 



Here /(q) is the Fourier transform of /y , while M{u!), 
reflecting the local dynamics, is given by the Weiss field 
Xo^ and the local spin susceptibility xioc as follows: 
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(6) 



Both Xo ^^'^ Xioc a-re entirely determined by the ef- 
fective impurity problem. The detailed diagrammatic 
derivations of Eqs. and (||) can be found in Ref. ||. 

A zero-temperature phase transition will occur in the 
Kondo lattice when xiQi^ = 0) becomes divergent for 
some momentum Q. In the extended DMFT formulation 
of the problem, M{uj) may be regular at the transition 
point [i.e., lmM{uj) ~ w], in which case the quantum 
phase transition will be of the usual spin-density-wave 
(SDW) type. However, should M{uj) instead have an 
anomalous frequency dependence at the transition, then 
the local critical behavior will also be anomalous. 

In order to distinguish between the two types of transi- 
tion, one must overcome the significant hurdle of solving 
a self-consistent Kondo problem with both a fermionic 
band and a vector bosonic bath. Progress in the stan- 
dard large--D DMFT was built on several decades' ac- 
cumulation of knowledge about the conventional Kondo 
and Anderson impurity models. Similar advances will be 
made within the extended DMFT only once the novel 
quantum impurity problem is understood. The next sec- 
tion outlines preliminary steps in this direction through 
analysis of the impurity problem without the imposition 
of self-consistency. 



III. QUANTUM PHASE TRANSITIONS IN 
GENERALIZED KONDO PROBLEMS 

Within the extended DMFT framework described in 
Section |l[ the nature of a quantum phase transition aris- 
ing in the Kondo lattice can be deduced from the fre- 
quency dependence of the quantity M{u}). Local Kondo 
physics will dominate if M{uj) is anomalous at the tran- 
sition. Since it is entirely determined by the effective im- 
purity problem [Eq. (||)], M{uj) can be anomalous only if 
the impurity problem has its own critical point with an 
anomalous local susceptibility. At the same time, self- 
consistency dictates that the spectral functions of the 
fermionic band and the vector bosonic bath must also 
be anomalous. Rather than determining the forms of 
these spectral functions using the extended DMFT self- 
consistency conditions, we have instead assumed certain 
unconventional forms from the outset. We show in this 
section that such forms can give rise to quantum critical 
points with local susceptibility anomalies. 



A. Kondo problem with an additional bosonic bath 

We first consider the Kondo impurity problem given 
m Eq. (|). Below some high-energy cutoff scale l/£,o, we 
assume that the conduction electron density of states is 
a (nonzero) constant. 




FIG. 2. Schematic renormalization-group flows for the 
Kondo problem with a additional bosonic bath described by 
an exponent A < 1. 
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and take the spectral function of the bosonic bath to be 
of the form 



q 



(8) 



In the self-consistent problem, A would be determined 
from the form of the local spin susceptibility Xioc- For 
a Fermi liquid, Imxioc('i^) ^ a; or equivalently, xioc(^) ^ 
as a result, the corresponding boson spectral func- 
tion would be Ohmic, i.e., A — 1. In a non- Fermi liquid 
for which xioc{t) ~ with a ^ 2, the boson spectral 
function would be non-Ohmic, i.e., ^ = a — 1 ^ 1. Here, 
we simply assume that A takes some fixed input value. 

The Kondo impurity problem described by Eqs. (|^), 
(^, and p^^— Jias recently been studied via a 1 — A 



expansion.OLL3 Upon progressive reduction of the high- 
energy cutoff, the couplings J and g entering Eq. (|^) 
satisfy the renormalization-group (RG) equations 

^ = («,/)- -if/, 



ding 
dln^ 



I- A 



1 



-{poJf-Kg^ 



(9) 



The RG flows for A < 1 are shown in Fig. ^. For 
small J and small g, there exists a separatrix specified 
by ^/Kg^ = (1 - A)exp[-(1 - A)/{poJ)]. The Kondo 
coupling J flows towards strong coupling for g < g^ 
and towards zero for g > g^. More generally, the 
flows are controlled by an unstable fixed point located 
at iVKg*,Po.n = iVT^, 1 - A). 

Within the 1 — A expansion, it is found the local spin 
correlation fupetion at the critical point has an anoma- 
lous exponentll3: x\oc (t) - l/tf^, where /3 = I - A. 
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B. Kondo problem with a pseudo-gap 

In the previous subsection, the fermionic spectral func- 
tion was assumed to be regular. It is possible, however, 
that critical fluctuations may produce a pseudogap in 
the effective single-particle density of states of the lattice 
model and, hence, in the density of states of the fermionic 
band of the impurity model. In order to isolate the effects 
of such a pseudogap, we neglect the bosonic bath in this 
subsection. Specifically, we study the SU(iV)-symmctric 
Kondo impurity Hamiltonian 



k,(T 



'COo 



(10) 



Here S is an impurity moment of spin degeneracy N, 
cocr (c = 1,2,..., N) annihilates an electron at the im- 
purity site, and r^^, (i = 1, 2, ... , A^^ — 1) is a generator 
of SU(A^). The conduction band density of states is as- 
sumed to take the power-law form 



Po|e| 




for |e| < 1, 
for lei > 1. 



fll) 



The problem described by EqSpdlO|) and ( |ll| ) was first 
studied by Withoff and Fradkin.EjThe local moment is 
quenched at low temperatures only if the Kondo coupling 
exceeds a threshold value Jc, where poJc ~ f- This makes 
the critical point at J = Jc a genuinely interacting prob- 
lem, in sharp contrast to the conventional (r = 0) Kondo 
model, for which Jc = 0. The strong-coupling (i.e., 
J > Jc) and weak-coupling (J < Jc) properties rof rthe 
power-law problem Jaave been studied extensively.Ey~t3 

In a recent work,E3 we studied the quantum critical 
behavior of the power-law Kondo model. In this context, 
it turns out to be crucial to distinguish between two dif- 
ferent spin susceptibilities: the thermodynamic quantity 
Ximp, which measures the response to a magnetic field 
that couples both to the impurity spin and to the con- 
duction electron spins, minus the equivalent response in 
the absence of the impurity; and the local susceptibility 
Xioc, the response to a magnetic field which couples only 
to the impurity spin. In the conventional Kondo problem, 
Xioc closely tracks Ximp as a function of temperature.Ell 
For a power law density of states, by caafa-ast, xioc and 
Ximp turn out to behave very differentlyEilcj as the Kondo 
coupling approaches Jc. Only xioc exhibits the scaling 
behavior characteristic of a continuous phase transition. 

This insight allows us to identify the form of the sin- 
gular part of the free energy, from which we derive scal- 
ing relations between the critical exponents in the quan- 
tum critical, weak-coupling, and strong-coupling regimes 
(see Fig. ||). Of particular interest is the exponent x 
which governs the temperature-dependence of the lo- 
cal susceptibility at the critical point via the relation 
Xioc{T,J = Jc) ~ 1/T^. As outlined below, x can be 
calculated in three separate limits. 



quantum 
*••., critical y'' 
weak \ /' strong 
coupling \ / coupling 
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FIG. 3. Schematic phase diagram of the Kondo problem 
with a power-law pseudogap. 

(1) For N ^ oo, the exponents on the strong-coupling 
side of the transition (J > Jc) can be calculated using 
slave-boson mean-field theory. Substituting these expo- 
nents into the scaling relations leads to the prediction 
that x{N — oo) — 1, independent of r. 

(2) When r is small, pqJc is also small, and x can be 
calculated directly via a perturbative expansion in po Jc. 
The result is 



x{poJc<t: 1) = 1- — (po^^c)^. 



(12) 



(3) For = 2, we have calculated the critical expo- 
nents using the numerical RG method. It is found that 
X varies continuously with r, and obeys Eq. ( p^ ) closely 
even for r as large as 0.45. 

The three independent methods for calculating x pro- 
duce results which are consistent with one another and 
point to the conclusion that for any finite iV, x varies 
continuously with r. 

In summary, we have studied two kinds of generalized 
Kondo impurity problems featuring unconventional host 
media. In both cases, there arises a quantum critical 
point with anomalous dynamics, such that the T = 
dynamical susceptibility varies as 



IniXioc(w) ~ (^^, 



(13) 



where B is less than 1 and depends on the form of either 
the boson spectral function or the fcrmion pseudogap. 



IV. QUANTUM PHASE TRANSITIONS IN 
HEAVY FERMIONS 

The results for generalized Kondo impurity prob- 
lems reviewed above, when combined with the extended 
DMFT formalism summarized in Section ||, suggest that 
it is indeed possible for a Kondo lattice system to display 
novel quantum critical behavior with anomalous local dy- 
namics. In this section we expand on this point. The ar- 
guments of this section, unlike those of the previous two 
sections, are necessarily suggestive in nature. 

We will find it useful to define an energy scale E^^^ 
below which M{uj) defined in Eq. (|^) is regular. Specifi- 
cally, 
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Im M (u;) = const, y— ^ for u < E^^, (14) 

where the constant is of order unity. 

For small values of the tuning parameter (5, the sys- 
tem lies well inside the paramagnetic region of Fig. ^. In 
such cases, the RKKY coupling Irkky is negligible, and 
the mean-field equations of the estended DMFT reduce 
to those of the standard DMFT .til The effective impu- 
rity problem is a Kondo model with a regular conduction 
electron band, for which both the local susceptibility x\oc 
and the Weiss field have the usual Fermi liquid form. 
Equation (|^) then implies that ImM(a;) ~ u>, i.e., that 
E^^^ is finite. 

Now let us consider larger values of (5, for which the 
RKKY interaction becomes more significant. Let us sup- 
pose that the extended DMFT for the Kondo lattice self- 
consistently generates an effective impurity problem with 
an unconventional fermionic or bosonic bath of the type 
discussed in the previous section. Varying the value of 5 
will change the couplings in the impurity problem. The 
results of Section III suggest that for some choice 5 — Jj^^., 
the impurity model will be tuned to a critical point char- 
acterized by a local susceptibility of the form given in 
Eq. dH). This in turn implies that at the critical point 
ImM(w) - cj^ with B < 1, and hence that £;j;^ ^ 0. 
We infer that E*^^ must vary with 5 roughly as shown in 
Fig. |. 

We can determine whether anomalous local physics 
plays an important role in the lattice quantum critical 
behavior by comparing 5^^^ defined above with 5c, the 
smallest value of 5 for which the momentum-dependent 
spin susceptibility x(Q,t<^ — 0) diverges at any wavevec- 
tor Q. 

If 8c is significantly smaller than S^oc^ then E^^^ re- 
mains finite (and acts as the effective Fermi energy) 
throughout the paramagnetic region right up to the criti- 
cal point. In this case, the magnetic phase transition rep- 
resents a Femri-surface instability induced by the RKKY 
interaction.E£l This is similar to the transition in more 
weakly intfimcting electron systems. The Hertz-Millis 
approachEEltJ should provide an appropriate basis for un- 
derstanding the quantum critical behavior. 

Alternatively, 5c may be equal (or very close) to Jj^j., in 
which case E*^^ vanishes (or nearly vanishes) at the crit- 
ical point. The quantum critical behavior in this case is 
dictated by the anomalous local dynamics. Equations (|^) 
and (|^) imply that will be anomalous over the en- 

tire Brillouin zone. 

The second type of critical point represents a continu- 
ous transition from a paramagnetic metal, having a Fermi 
surface of volume Nf -\- Nc, to a magnetic metal, with a 
Fermi surface of volume Nc- Here Nf and Nc denote the 
number of /-electron spins and conduction electrons, re- 
spectively. One can think of the anomalous local dynam- 
ics as capturing the fluctuations between the two Fe™i 
surfaces. (This last picture has been noted previously.c3) 

The discussion to this point has focused on the inter- 
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FIG. 4. The local energy scale defined in Eq. (^_^ as a 
function of the tuning parameter. 

mediate quantity, M{uj). However, the measurable quan- 
tity in a Kondo lattice system is x{^.^ From Eq. (||) it 
is apparent that the low-frequency behavior of ImAf(a;) 
is the same as that of Im xioc , where 

generic q 

Here "generic q" restricts the summation to wavevectors 
that are not too closecJ to the ordering wavevector Q. 

The preceding discussion suggests an operational crite- 
rion for distinguishing the two types of critical behavior. 
Let us redefine the energy scale E*^^ so that 

Imxioc^t^ for ^ < -EToe- (16) 

Then, if E^^^ = at the quantum critical point, the phase 
transition is driven by local physics. Otherwise, the tran- 
sition is of the usual SDW type. 

Neuteoij^attering experiments have recently been car- 
ried outlla'E2lon CeCug-xAua; at the critical concentration 
Xc = 0.1. At wavevectors close to the ordering wavevec- 
tors, the dynamical spin susceptibility appears to have 
an anomalous frequency dependence, with an exponent 
B K, 0.8. Furthermore, the same exponent also describes 
the q = susceptibility, suggestingfij that the anomalous 
energy exponent occurs over a large region of the Bril- 
louin zone. The latter observation, according to Eqs. ( |T5| ) 
and (|l6|), implies that E^^^ = at the critical point. 
Thus, the quantum critical behavior of CeCug-xAua; ap- 
pears to fall in the category of quantum phase transitions 
driven by local physics. 

V. SUMMARY AND OUTLOOK 

To summarize, we have reviewed an extended DMFT 
which maps the Kondo lattice problem onto a gener- 
alized Kondo impurity problem, supplemented by self- 
consistency conditions. The Kondo and RKKY physics 
of the lattice model are both manifested in the effective 
impurity problem: The former is described by the cou- 
pling of the local moment to a self-consistent conduction 
electron band and the latter by the coupling of the local 
moment to a self-consistent bosonic bath. 
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Two kinds of quantum phase transitions are argued to 
occur in Kondo lattice systems. One is the usual SDW 
type, whereas the other is governed by anomalous local 
Kondo dynamics. For quantum critical behavior of the 
second type to occur, it is necessary that the generalized 
Kondo impurity arising in the extended DMFT displays 
a critical point at which the local susceptibility has an 
anomalous exponent. Such critical points occur in two 
related Kondo impurity models, one describing a local 
moment coupled to an additional bosonic bath, the other 
featuring a pseudogap in the conduction band. 

Finally, we have introduced an operational criterion 
that can be used to identify quantum critical behavior 
driven by local physics. It is based on the vanishing of 
a local energy scale [E;^^^ defined in Eqs. and ^] 
at the critical point. The results of neutron scattering 
experimentalj suggest that this criterion is satisfied in 

CeCu6-2;AU:E. 

A number of significant issues remain to be addressed. 
Most importantly, the self-consistent mean-field equa- 
tions need to be analyzed more completely. Once that 
has been achieved, it should be possible to study the 
effect of disorder on this novel type of quantum phase 
transition.Eil Finally, the effect of valence fluctuations can 
also be addressed within this framework. 

We would like to thank P. Coleman, G. Kotliar, A. 
Rosch, S. Sachdev, A. Sengupta, and C. M. Varma for 
useful discussions. This work has been supported in part 
by NSF Grant No. DMR-9712626, Research Corpora- 
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